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Abstract 

The PDFs for energy dissipation rates created in a high resolution from 4096"^ DNS for fuUy devel- 
oped turbulence are analyzed in a high precision with the PDF derived within the formula of multifractal 
probability density function theory (MPDFT) . MPDFT is a statistical mechanical ensemble theory con- 
structed in order to analyze intermittent phenomena through the experimental PDFs with fat-tail. By 
making use of the obtained w-PDFs created from the whole of the DNS region, analyzed for the first time 
are the two partial PDFs, i.e., the max-PDF and the min-PDF which are, respectively, taken out from 
the partial DNS regions of the size 512"^ with maximum and minimum onstropy. The main information 
for the partial PDFs are the following. One can find a w-PDF whose tail part can adjust the slope of 
the tail-part of a max-PDF with appropriate magnification factor. The value of the point at which the 
w-PDF nmltiplicd by the magnification factor starts to overlap the tail part of the max-PDF coincides 
with the value of the connection point for the theoretical w-PDF. The center part of the min-PDFs can 
be adjusted quite accurately by the scaled w-PDFs with a common scale factor. 



1 Introduction 



There are the keystone works [Tl[ll3ll[51ini[71[Hl[i[ini[ni[Ill[IllIl[IIllII| providing the multifractal 
aspects for fuUy developed turbulence. Most of the works [DiaElllElElHElIinilll] deal with the 
scaling property of the system, e.g., comparison of the scaling exponents of velocity structure function. 
Only a few [51 [131 [IH [13 [I^ analyze the probability density functions (PDFs) for physical quantities 
representing intermittent character. Among these researches, multifractal probability density function theory 
(MPDFT) [HI O [H [13 dl] provides the most precise analysis of fat-tail PDFs, which is a statistical 
mechanical ensemble theory constructed by the authors (T.A. and N.A.) in order to analyze intermittent 
phenomena through fat-tail PDFs extracted out from the data by experiments or numerical simulations. 

In order to extract the intermittent character of the fully developed turbulence, it is necessary to have 
information of hierarchical structure of the system, which is realized by producing a series of PDFs for 
responsible singular quantities with different lengths 

l^^i^S-^, <5>1 (n = 0,1,2, •••) (1) 

that characterize the regions in which the physical quantities are coarse-grained. The value for 5 is chosen 
freely by observers. NETFD is constructed such a way that the choice of 6 should not affect the theo- 
retical estimation of the values for the fundamental quantities characterizing the turbulent system under 
consideration. 

Within MPDFT, it is assumed that there are two contributions to form a fat-tail PDF, i.e., one is a 
coherent contribution coming from a coherent turbulent motion and the other is a incoherent contribution 
from the dissipation term in the N-S equation which violates the invariance under the scale transformation 
given at the beginning of the next section. It is also assumed that the fat-tail PDF can be divided into two 
parts, i.e., the tail part and the center part. As the tail part is responsible for the intermittent character 
of the system, the coherent contribution may dominate over the incoherent contribution at the tail part. 
On the other hand, the center part consists of both coherent and incoherent contributions. With the above 
consideration, we neglect within A&A model of MPDFT the incoherent contribution to the tail part in 
deriving the theoretical PDF formula, and conjecture that the coherent contribution to the tail part can be 
determined solely by the PDF for a provided in the next section. Through the conjecture, we are able to 
separate theoretically the incoherent contribution from the coherent one at the center part. 

In this paper, we analyze in a high accuracy the PDFs for energy dissipation rates created in a high 
resolution from the snapshot raw data of 4096^ DNS for fully developed turbulence conducted by the authors 
(Y.K. and T.I.) [19^. By making use of the result obtained by the precise analyses on the PDFs obtained from 
the whole of the DNS region, we analyze for the first time those PDFs extracted from the partial DNS regions 
of the size 512'^ obtained by cutting the whole of 4096"^ DNS region into 512 pieces. In Section [2] a brief 
introduction of A&A model is given with its concept and basics. In Section[3| the formula for the theoretical 
PDF for energy dissipation rates is derived within A&A model of MPDFT. The two distinct divisions of the 
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PDF are introduced, which are important in the following analyses. In Section |4j the high-resolution PDFs 
created from the whole of 4096'^ DNS region are analyzed in a high precision with the derived theoretical 
PDF. In Section [sj the partial PDFs created from the 512^ partial DNS regions with maximum enstrophy 
and minimum enstrophy are studied with the help of the PDFs extracted from the whole of 4096"^ DNS 
region. Conclusion is provided in Section [6j 

2 A&A model 

It is known that the Navier-Stokes (N-S) equation du/dt + {u ■ V)m = —Vp + vV^u for an incompressible 
fluid (V • M = 0) is invariant under the scale transformation x Xx accompanied by the rescaling in velocity 
field, time and pressure, i.e., u — > A'^/'^u, t — > X^^^^^t and p — >■ X^^^^p with an arbitrary real number a, 
when V = Q. Here, v is the kinematic viscosity; u is the velocity field; p is the pressure of fluid per mass 
density. It is assumed that for high Reynolds number {v — >■ 0) the singularities distribute themselves in a 
multifractal way in real physical space, and that they are the origin of a coherent turbulent motion providing 
the fat-tail part of PDFs [2 . In treating an actual turbulent system, the value of v is fixed to a finite non- 
zero value unique to the material of fluid prepared for an experiment. Therefore, for the study of fully 
developed turbulence, we have to look for the characteristics of a coherent turbulent motion surrounded by 
an incoherent fluctuating motion. The latter motion is due to the dissipation term vXI^u in the N-S equation 
which is the term violating the invariance under the scale transformation. MPDFT is a statistical mechanical 
ensemble theory for analyzing both the coherent turbulent motion and the incoherent fluctuating motion. 
The dissipation term can become effective depending on the region under consideration since the spots of 
the region where the invariance is broken locate here and there in fully developed turbulence, and provide 
us with a speciflc character of fluctuation around the coherent motion. 
Let us consider the energy dissipation rate e„ that is related to a by 

Sn = . (2) 

Note that lim„_j.oo £„ ^ oo for a < 1. The degree of singularity for energy dissipation rate is specified by 
the singularity exponent a which is considered to be a stochastic variable whose PDF, P(")(q;), is given 
by Uni [n H 113 III EZl [H] the Renyi [ID] or HCT ^23 Hg type PDF 

p(»)(«)oc[l-(a~ao)V(^«)T^''"'^ (3) 

with Aa = ^2X/(1 - q) In (5. This is the MaxEnt PDF of the Renyi entropy or the HCT entropy The 
domain of a is amin < a < amax with amin and amax being given by amin/max = Q^o T ^Q?- <Z is the entropy 
index. The multifractal spectrum introduced through P("^(q;) cx (^n/^o)^"''^''"-' is given for n 1 by 

/(a) = 1 + {In [1 - (a - ao) V(^«)'] } /(I - In 5. (4) 

^ The function ([3| is the MaxEnt PDF derived from the Renyi entropy or from the HCT entropy with two constraints, one 
is the normaUzation condition and the other is a fixed g- variance |22l . 
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The three parameters ao = (a), X and q are determined by the following three conditions where the 
average (• • •) is taken with P*^"''(a). One is the energy conservation law (e„) = e. Another condition is 
the definition of the intermittency exponent /i, i.e., ((e„/e)^) — (i?„/^o) The third condition is the new 
scaling relation 

ln2/(l-(7)In(5 = l/a_-l/a+ (5) 

with a± being the solution of f{a±) = 0, i.e., a± = ao ± {2hXY/^ with b = (1 - c-(i-«)^"'')/(l - q)hi5. 
The scaling relation ^ is solved to give (26X)i/2 = -(1 - q) \og^ 5 + {al + [(1 - q) log2 (5]^}^/^. The new 
scaling relation is a generalization of the one introduced by Tsallis and others [231 HI] to which ^ reduces 
when 5 = 2. This generalization was born out of A&A model within the theoretical framework of MPDFT 
itself. Note that the A&A model is a one-parameter model depending only on /i. 

The parameter q is determined, altogether with ao and A", as a function of ^ only in the combination 
(1 — g) In (5. It is quite reasonable in the following reason. The value of the magnification 5 is determined 
arbitrarily by observers, therefore its value should not affect the values of physical quantities as long as one 
studies a single turbulent system. The difference in 5 is absorbed into the entropy index q, therefore changing 
the zooming rate 5 may result in picking up the different hierarchy, containing the entropy specified by the 
index q, out of multifractal structure of turbulence. 

3 PDFs of energy dissipation rates with A&A model 

Let us create the probability n^"^(e„)(ie„ to find the physical quantity e„ taking a value in the domain 
£n ~ £n + dem whosc normalization is specified by JJ^ denn*-"^ (£„) = 1. We assume that the PDF, n(")(e„), 
can be divided into two parts as 

n(")(e„) = n(")(£„) + An(")(£„). (6) 

The first term is the part representing a coherent turbulent motion described in the limit — > 0. The 
coherent contribution is assumed to be given by 

n(")(|e„|)d|e„|=n(")p(")(a)da (7) 

with the variable translation ^ between e„ and a. On the other hand, the second term An'^"^(e„) in (joj) 
represents an incoherent contribution from the dissipation term in the N-S equation. The dissipation term 
violates the invariance under the scale transformation, and therefore the incoherent contribution has not 
been taken into account in the consideration given in the previous section for P'^"^(a), i.e., the effect of the 
finiteness of the dissipation term is not included in ng"^(e„). Each term is composed of the product of two 
PDFs, i.e., 1) the PDF that determines from which the contribution is originated out of two independent 
origins, i.e., the coherent origin or the incoherent origin, and 2) the conditional PDF to find a value £„ in 
the domain e„ ^ e„ + rf£„ for each origin. 
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We divide the PDF n("^(e„) into two parts another way, i.e., 

n(")(£„) = n(;^)(£„) + n[r^(e„). (8) 

The tail part n|"''(e„) for |£„| > e* (equivalently, < « < ct*) represents mainly the contribution 

of the singularities due to the coherent turbulent motion, whereas the center part n^^"|.(e„) for |e„| < e* 
(equivalently, a,nin > < «*) does the contributions from both coherent and incoherent motions. Here, e* 
is the connection point of the two PDFs, 11^" •*(£„) and n|"''(£„), which is related to a* through It is 
reasonable to assume that the origin of intermittent rare events is attributed to the first singular term in ^ , 
and that the contribution from the second term Z\n^"^(e„) to the events is negligible. We assume therefore 
that, for the tail part PDF n[f^(£„), 

one can neglect completely the contribution from the second correction 
term in (|6]). Under this assumption, we put 

nlr'(£„)d£„ = n(")p(")(a)da (9) 



for £„ > £;. 

Let us introduce, here, the PDF Il(")(^„), through the relation H'-"-' (Cn)c?Cn = n(")(£„)(i£„, for the new 
variable f„ = £„/ {{ef-^))l^'^ normahzed by the standard deviation. Here, (((£,i)^))c = (((sn)^)) — ((^n))^ ^i^"^ the 
average ((• • •)) is taken with n'^"-'(£„). The normalized variable ^„ is related to a by ^„ = fnl^n/^o)"^^"^^ 
with Ca™ = 1 - r(m) and e„ = [{7^"^ " [tJ"' + (1 - 7^"V3(C/4)'^-in(4./4)-'^« + (1 - T^'^Ve]"'/" where 
a3m = (|/"(«o)//"(a™)|)^/^7i"' = (^„/4)'"/„°°d£„(£„)"An(")(£„)andT(q-) = l-aoq+2Xqy{l+C'/') + 
[1 - log2(l + C^/^)](ln2)/(l -q)\nS with Q = 1 + 2^2^(1 - q) \n5. The mass exponent T{q) is related to 
f{a) by the Legendre transformation [31[S] /(a) — aq + T{q) with q — df(a)/da and a — —dT{q)/dq. The 
variables q and a are related through q — [2/(1 — q) In J] [(a — aQ)/{a — amax)(<^ ~ ckmin)]- The gth moment 
of the energy dissipation rate is given by ((£„/e)«) cx (4»/4)"^''''+^"'^ I3[I71[IHj. With the assumption 
we have, for Cn < \U < C (equivalently, a,„in < a < a*) with C"' = 6.(^«/4)"»"-'^«/^ 

n[r^ (e„) = n(") (^„/4)'"^^"^ U^n (10) 

where H'") = nf,"^f„(|/"(ao)|/27r| ln(C/4)|)i/2 with H^"^ = 1 - 7^"^ Here, ^ and a* are related with each 

other hy C = Caila/ior'-^"^^- 

It may be appropriate to perform the connection by means of tl'-^\^„) and n|r^(^„). The tail and 
the center parts of PDFs are connected at ^* under the conditions that they have the common value, i.e., 
n^rHC) = Ur\Cn) and the common log-slope, i.e., dlnll[^\^r.) / dUir.=C„ - dlnni?)(e„)/d^„le„=e. . The 
value ^* (therefore, £* ) is determined for each PDF as an adjusting parameter in the analysis of PDFs, 
and provides us with an important information of the border line beyond which the contribution from the 
incoherent motion ceases. 
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As there is no theory at present to produce the formula of fli^\^,n) for < ^* (equivalently, a* < a), 
we introduce a trial PDF which, after the connection procedure, is given by 

fi^'HU = n(") (ejioy-f^'^'^ e-^9iu)~9icj] (11) 

with a trial function of the Tsallis-type e-9(C.O = (^^/^* {1 - (1 - q') [e + /'(a*)] [(Cn/C)"' - 1] /wf^^^^' 
and n'"^ = ng"' ■\/|/"(Q!o)|/27r| ln(4i/^o)l/Cn- The parameter w is adjusted by the property of the exper- 
imental PDFs around the peak point; q' is the entropy index different from g in ([s]); is determined by 
the property of PDF near ^„ = 0. The contribution to nCT''(^„) comes both from coherent and incoherent 
motions. 

Thanks to the assumption that, for the tail part ^* < ^„ < S,™'^^, one can neglect the contribution from 
An(")(a;„), we obtain the formula to calculate 7™ ■* in the expression 7^ "* = (-^m ■* — Lm^)/ (1 + K^'^ — £q"'') 
with 

^i"^ = (^«/4)'^^("'^+'""Vl/"(«o)|/27r|ln(C/4) I C dz z™e-[s(«'*-)-»«*.)l , (12) 
Lt^ = (^n/4)Vl/"(«o)||ln(£„/4)|/27r j""^^^ da (13) 

J a* 

4 Analysis of PDFs taken from the whole of 4096^ DNS region 

PDFs of energy dissipation rates studying in this section are extracted by cooking the snapshot data taken 
from the whole of 4096'^ DNS region. The DNS [H] is specified by the Taylor micro-scale Reynolds number 
R\ ~ 1132 and the integral length L/77 = 2.130 x 10^ with the Kolmogorov length scale 77 = 5.12 x 10^"*. 
The inertial range is estimated to be in the region 126 < r/rj < 448. 

The PDFs are analyzed in Fig. [l] for three different magnifications (a), (b) 5 = 2^/^, (c), (d) 5 — 2^/^ 
and (e), (f) (5 = 2. The vertical axes of (a), (c), (e) in Fig. [I]are given in log scale which are good to see the 
tail parts of PDFs, whereas those of (b), (d), (f) in the figure are given in linear scale which are appropriate 
to observe the center parts of PDFs. For better visibility, PDFs in each figure are shifted, successively, with 
respect to the values of r/r] by appropriate amounts along the vertical axis, i.e., —1 in (a), —0.1 in (b), —2 
in (c), —0.2 in (d), —4 in (e) and —0.4 in (f). The PDFs with smaller values of r/77 are placed at upper 
parts in each figure and therefore PDFs with larger values of r /r] at lower parts. Open circles in the figures 
are the PDFs constructed with the help of the fiuid velocity data taken from all the points in the whole of 
4096^ DNS region. Solid lines represent the theoretical PDFs. 

In 18J, the authors (N.A. and T.A.) performed, by means of the theoretical PDF, the analysis of the 
PDFs furnished by Kaneda and Ishihara [19 which had been created by setting the width of bins to be 0.2 
along the ^„ axis. The resolution of the PDFs is good enough to analyze the tail part but is not enough to 
analyze the center part (see Fig. 1, Fig. 2 and Fig. 3 in [18 ), which results in the difficulty in drawing precise 
theoretical curves for the center-part PDFs. This ambiguity causes the scattering of parameters, especially. 
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Figure 1: PDFs of energy dissipation rates (a), (b) for S = 2^/'*, (c), (d) for S = 2^/^ and (e), (f) for 6 = 2. 
Open circles represent the PDFs extracted out from the whole of 4096'^ DNS region. For better visibility, the 
PDFs in each figure are shifted along the vertical axes by the amount —1 in (a), —0.1 in (b), —2 in (c), —0.2 
in (d), —4 in (e) and —0.4 in (f). They are displayed according to the values of r/rj (= in/v) successively 
(see Table [T]) in the sequence that the PDF with the smallest value of r/rj is placed at the top and the one 
with the largest value at the bottom for each 6. Solid lines represent the curves given by theoretical PDFs 
of A&A model. It turns out that n = 0.345 ((1 - q)\nS = 0.299, ao = 1.20, X = 0.411) for every S. Other 
parameters are listed in Tablejl] The vertical axes in (a), (c), (e) are of log-scale, whereas those in (b), (d), 
(f ) are of linear-scale. 
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Table 1: The values of parameters in the theoretical PDF which are obtained in the course of the analyses 
of PDFs extracted out from the whole of 4096^ DNS region for the cases S = 21/-*, S = 2^/2 and (5 = 2. 
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27.5 


29.0 


5.03 


1.06 


0.480 


2.90 


0.220 


2.84 


14.0 


4.85 


1.05 


0.450 


2.90 


0.210 


2.03 


8.00 


5.55 


1.05 


0.4.50 


2.90 


0.220 


2.28 


32.8 


28.0 


4.85 


1.06 


0.490 


3.10 


0.210 


2.83 






























38.9 


26.5 


4.59 


1.06 


0.500 


3.00 


0.200 


2.39 


13.5 


4.68 


1.06 


0.510 


3.00 


0.210 


2.28 
















46.3 


26.0 


4.51 


1.06 


0.320 


3.30 


0.200 


2.32 






























55.1 


25.0 


4.33 


1.08 


0.600 


3.00 


0.200 


2.16 


12.5 


4.33 


1.08 


0.600 


3.00 


0.200 


2.17 


6.30 


4.37 


1.07 


0.600 


3.00 


0.200 


2.09 


65.6 


24.0 


4.16 


1.07 


0.600 


3.40 


0.200 


1.97 






























78.0 


23.0 


3.99 


1.08 


0.650 


3.30 


0.200 


1.68 


11.5 


3.99 


1.08 


0.650 


3.30 


0.200 


1.68 
















92.7 


22.0 


3.81 


1.08 


0.650 


3.70 


0.200 


1.49 






























110 


21.0 


3.64 


1.08 


0.700 


3.80 


0.200 


1.34 


10.7 


3.71 


1.08 


0.700 


3.80 


0.200 


1.34 


5.30 


3.67 


1.08 


0.720 


3.80 


0.200 


1.38 


131 


20.0 


3.47 


1.08 


0.730 


3.80 


0.200 


1.09 






























156 


19.0 


3.29 


1.09 


0.770 


3.85 


0.190 


0.997 


9.5 


3.29 


1.08 


0.750 


3.90 


0.200 


0.920 
















186 


18.0 


3.12 


1.09 


0.800 


3.90 


0.200 


0.808 






























221 


17.0 


2.95 


1.09 


0.850 


4.10 


0.190 


0.733 


8.50 


2.95 


1.09 


0.850 


4.30 


0.200 


0.700 


4.30 


2.98 


1.08 


0.780 


4.60 


0.180 


0.7.53 


264 


16.0 


2.77 


1.10 


0.900 


4.20 


0.180 


0.681 






























314 


15.0 


2.60 


1.10 


0.920 


4.50 


0.180 


0.559 


7.50 


2.60 


1.09 


0.900 


4.. 50 


0.180 


0.570 
















374 


14.0 


2.43 


1.09 


0.950 


4.70 


0.180 


0.473 






























442 


13.0 


2.25 


1.10 


1.02 


4.85 


0.160 


0.448 


6.50 


2.25 


1.10 


1.00 


5.00 


0.160 


0.458 


3.30 


2.29 


1.10 


1.05 


4.70 


0.170 


0.438 



in the r ///-dependence of 9 (see Fig. 7 in |18j). Note that r {— £„) in the present paper corresponds to 2r 
(= C) in dlllH]. 

In order to raise the resolution of PDFs, in the present paper, we created PDFs by cooking the row data 
of 4096"^ DNS turbulence. In creating the tail-part PDFs in Fig. [l] we set the width of bins to be 5 x 10~^ 
along the ^„ axis, while, in creating the center-part PDFs, we set the width of bins to be 5 x 10^"^. Note that 
in drawing the tail-part (center-part) PDFs, not all the bins but every 20 (10) bins are plotted for better 
visibility. We discarded the bins containing the number of data points less than 10^^ % of the mean number 
of data points per bin on average. For example, the bins containing less than 1.37 x 10^ (9.73 x 10^) data 
points are discarded for r/rj — 27.5 {r/ri = 443). Note that there are 4096^ = 6.87 x 10^° data points as a 
whole. 

The parameters necessary for the theoretical PDF within A&A model, i.e., those for the tail part of 



the PDF (10 1 and for the center part of the PDF (111, are obtained through the analysis of the high 
resolution PDFs created from 4096'^ DNS. It turns out that the tail-part PDFs for the turbulent system 
under consideration are characterized with the value of the intermittency exponent /i = 0.345. Then, the 
parameters necessary for the PDFs are determined as [1 — q)h\5 ~ 0.299, — 1-20 and X ~ 0.411, which 
are independent of 5. It may be worthwhile to note here that the entropy index becomes q = —0.728 for 
5 ^ 2^/^ (= 1.19), q = 0.136 for 5 = 2^/'^ (= 1.41) and q = 0.568 for 5 = 2. The parameters necessary for 
the center-part PDFs are listed in Table [T] for each 5. Note that the value of /i extracted from the PDFs 
with high resolution in the present paper and those with low resolution in |18| turns out to be the same, and 
there is no difference observed between the tail-part PDFs of two resolutions. 

Fig. [2] and Fig. [3] give the dependence of n, q' , 9, Inw and ln£„ on r/77 (= in/r]) for S = 2}!^ (A), 
for (5 = 2^/^ (x) and for (5 = 2 (Q) extracted from the series of PDFs (see Fig. [I]). Here, h is defined by 
n = nln(5 with which introduced in ([l]) reduces to £„ = €oe~". The solid lines in Fig. [2] (a), (b), (c), (d) 
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Figure 2: The r/rj {= in/v) dependence of the parameters (a) n, (b) q' , (c) 9 and (d) w for the PDFs of energy 
dissipation rates, extracted from the whole of 4096'^ DNS region. In each figure, the resuhs for S = 2^/^ (A), 
5 = 2^/^ ( X ) and 5 = 2 (Q) are displayed altogether. The line in each figure is obtained by the method of least 
squares using the data points for every S. The obtained formula is (a) n = —2.30 ^ogiQ{r/r]) + 8.37, (b) q' = 
3.71x10-2 logio(r/7;) + 1.00, (c) 9 = 1.54 logio(r/?7) + 0.585 and (d) logio w = 0.280 logio(r/77)+logio 0.186. 
The inertial range is the region between the vertical dash-dotted lines. 
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Figure 3: The r/rj {= in/rj) dependence of the connection point £* for PDFs of energy dissipation rates, 
extracted from the whole of 4096^ DNS region. The results for 5 = 2^/^ (A), 5 = 2^''^ (x) and (5 = 2 (Q) are 
displayed altogether. The hne log^Q e* — —0.727 logio(»'/'7) + logio 37.8 is obtained by the method of least 
squares using the data points for every 5. The inertial range is the region between the vertical dash-dotted 
lines. 



and in Fig. [3] represent the empirical formulae given, respectively, by 

n = -0.998 ln(r/77) + 8.37, 

g' = 1.61 X 10^2 ln(r/?7) + 1.00, 

d = 0.668 ln(r/77) + 0.585, 

\nw = 0.280 ln(r/77) - 1.68, 

ln4 = -0.727 ln(r/77) + 3.63, 



(14) 
(15) 
(16) 
(17) 
(18) 



which are obtained by the method of least squares using in each figure all the data points for 6 — 2^/^, 2^/^ 
and 2 altogether. The results given in Fig. [2] and Fig. |3] prove the correctness of the assumption that the 
fundamental quantities of turbulence are independent of 5. Note that, in the captions of these figures, we 
chose the base of logarithmic function in the empirical formulae to be 10 which corresponds to the axes of 
figures. 

With the PDFs of high resolution, especially, for the center part given in Fig. [T] we succeeded to get 



the correct empirical formulae (14), (15), (16), (17) and (18). The independence of h from 5 ensures the 



uniqueness of the PDF of a for any value of 5 when the intermittency exponent /i has been settled. There 
appears big difference in the parameters q', 9 and w responsible for the center-part PDFs compared with 
the results in [18 obtained by PDFs with low resolution. It turns out that q' ~ 1.08 but depending slightly 
on r/rj with positive slope for the high resolution (see Fig. [2](b)), whereas it had negative slope for the low 
resolution (see Fig. 6 in |18j). The scattering of 9 in the dependence oir/rj has been reduced (compare Fig. [2] 
(c) with Fig. 7 in [H]). It is found that Inw depends lineally on ln(r/77) which is observed also in the analysis 
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Figure 4: Comparison of the max-PDF (□), the min-PDF (A) and the w-PDF (Q) for r/77 = 55.1. The 
horizontal axes represent the energy dissipation rates £„, and the vertical axes do the PDFs on (a) log and 
(b) linear scale. 



of PDFs for energy dissipation rates extracted from the turbulence in a wind tunnel |25j , while w was almost 
constant in the analyses of [18j with less resolution. The connection point a* is adjusted in order for the 
best fit of the PDF around the region between the peak and the connection point. Note that the region 
a < a* {a > a*) corresponds to the tail (center) part of the PDF. It is revealed that the value a* satisfies 
a* ~ 0.2 for all the data points with different values of r/rj (see Table [T]), which proves the assumption that 
the center part nCT^(£„) is constituted by two contributions, one from the coherent contribution ng"'(£„) 
and the other from the incoherent contribution An'^"^(e„), and that almost all the contribution to the tail 
part n|"''(e„) comes from the coherent motion of turbulence. Remember that the energy dissipation rate 
becomes singular for a < 1. 



The comparison of the extracted formula (14) for h with the theoretical relation h = —\n{r/r])+ln{io/ri) 
provides us with the estimation io/rj — 4.31 x 10'^. Since the smallest grid spacing is 3r/ 19J, lo/irj = 1.44 x 10'^ 
provides us with the number of grids corresponding to £o- Note that (q/t] is about 2 times larger than the 
integral length L/rj oi the system [12]. It should be noted here that we observed £o/ri ~ L/rj for the case of 
experimental turbulence in a wind tunnel [25j . 

5 Analysis of PDFs taken from partial DNS regions 

We are analyzing in this section the PDFs of energy dissipation rates created from the snapshot data in 
partial regions of the size 512'^ which are obtained by cutting the whole of 4096^ DNS region into 512 pieces. 
We will refer to the PDF created from the partial region as p-PDF in the following. Among 512 partial 
DNS regions, we select in this paper two DNS regions, i.e., one has a maximum cnstrophy and the other a 
minimum enstrophy, and study the series of p-PDFs obtained in these regions. We call the p-PDF created 
from the partial DNS region with maximum (minimum) enstrophy max-PDF (min-PDF). 

The max-PDF (□) and the min-PDF (A) for the case of r/rj — 55.1 are displayed as functions of e„ in 
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Figure 5: Analyses of (a) the max-PDFs (□) and (b) the min-PDFs (A) of energy dissipation rates for 
,5 = 2^/"*. For better visibihty, the PDFs in each figure are shifted along the vertical axes by the amount — 1 
in (a), -4 in (b) according to the values r/rj = 27.5, 32.8, 38.9, 46.3, 55.1, 65.6, 78.0, 92.7, 110, 131, 156, 
186 and 221 for (a), and r/r/ = 27.5, 32.8, 38.9, 46.3, 55.1, 65.6, 78.0, 92.7 and 110 for (b), successively, from 
the top to the bottom. Solid lines are the w-PDFs drawn in Fig. [2] after appro pria te cooking procedures. 



i.e., (a) multiplying some factor X to w-PDF, (b) scaling w-PDF as n(")(£„) in (n9|, whose recipe is given 
in Fig. [6] (see the main text for detail). 



Fig. |4]with the vertical axes in (a) log scale and (b) linear scale. In the figure, we put the PDF (Q) created 
from the whole of 4096'^ DNS region, which we call w-PDF in the following, for the same value of r/rj as 
a reference. We observe by comparing with the w-PDF that the proportion of the probability density for 
the max-PDF (the min-PDF) is shifted to the tail-part PDF (the center-par PDF). It is reasonable in the 
sense that since the vortexes are distributed dense (sparse) in the partial region with maximum (minimum) 
enstrophy, the proportion of the intermittent coherent (the fluctuating incoherent) fluid motion should be 
large. 

The max-PDFs (□) and the min-PDFs (A) of energy dissipation rates for 5 — 2^/^ are given, respectively, 
in Fig. [5] (a) and (b) . For better visibility, the PDFs in each figure are shifted along the vertical axes by 
the amount -1 in (a), -4 in (b) according to the values r/-q = 27.5, 32.8, 38.9, 46.3, 55.1, 65.6, 78.0, 92.7, 
110, 131, 156, 186 and 221 for (a), and r/rj = 27.5, 32.8, 38.9, 46.3, 55.1, 65.6, 78.0, 92.7 and 110 for 
(b), successively, from the top to the bottom. Solid lines are the w-PDFs drawn in Fig. [2] (a), (b) after 
appropriate cooking procedures whose recipe is listed in Fig. [6| From Fig. [5] (a), it was revealed that the 
tail-part of max-PDFs can be adjusted with the slope of the tail-part of w-PDF with a specific value of r/rj 
magnified by some factor x, and that the value of e\ at which a max-PDF (□) and a solid line representing 
a w-PDF in Fig. [5] (a) start to overlap is quite close to the connection point e* associated to the w-PDF 
(see Fig. |6|. From Fig. [5] (b), it was found that the center-part of min-PDFs can be adjusted by the scaled 
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Figure 6: Recipe to analyze the max-PDFs and the min-PDFs by means of the w-PDFs. The arrows indicate 
that which w-PDF with a specific r/77 at the origin of an arrow can adjust the max-PDF or the min PDF at 
the point of the arrow, respectively, by the magnification factor x and by the scaling factor a. is the point 
from which a max-PDF (□) and a solid line representing a w-PDF in Fig. [s] (a) multiplied by the factor x 
start to merge. The values of e* associated with w-PDF for S = 2^/'* are listed here again for convenience. 



PDF, n(")(e„), which is introduced through 

n("n£n) den = n("H£n) dSn, £„ = £„/a (19) 

with an appropriate scaling factor cr (see Fig. [6]) . 

In Fig. [6j we listed the values of the magnification x which is necessary to analyze the max-PDFs in Fig. [5] 
(a), and of the scaling factor a which is necessary for the scaled PDF Il(")(£„) to analyze the min-PDFs 
in Fig. [5](b). The arrows indicate that which w-PDF with a specific value r/i] at the origin of an arrow 
can adjust the max-PDF or the min-PDF at the point of the arrow. In between the arrows, there inserted 
the magnification factor x for the case of max-PDF or the scale factor a for the case of min-PDF (see also 
Fig.§. 

From Fig. |6j we notice at least two remarkable outcomes. One is the fact that the value of ej^ associated 
with max-PDF is almost equal to the values of the connection point e* , which may indicates that the division 
of the tail part and the center part of PDF within A&A model is reasonable. The other is the fact that the 
scaling factor cr has almost a common value around 2.19, which may indicate that there exists a beautiful 
statistics associated with the incoherent fluctuation part of turbulent motion. The latter may be related 



to the beautiful scaling behaviors given by the empirical formulae (15), (16) and (171 for the parameters 
associated with the center-part PDFs. 

6 Conclusion 

We investigated quite accurately the PDFs of energy dissipation rates created in a high resolution by cooking 
the snapshot data taken from the whole of 4096^ DNS region by the theoretical formula for PDF derived 
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within A&A model of MPDFT whose contents are compactly given in the present paper. Analyzing the 
obtained high-resolution w-PDFs, we derived the empirical formulae of the parameters consisting, especially, 
the central part of the theoretical PDF in their most precise forms for the first time. It was revealed that 
n, q' , 6, \nw and Ine* (therefore, a*) are independent of S thanks to the new scaling relation ([5|, and that 
they show scaling behaviors extending to the regions with smaller r/rj values from the inertial range. By 
making use of the w-PDFs, we also succeeded to extract the some attractive informations contained in the 
max-PDFs and the min-PDFs (see Fig|6]): 1) We can find a w-PDF whose tail part can adjust the slope of 
the tail-part of a max-PDF with appropriate magnification factor %. 2) The value of ejj at which the w-PDF 
multiplied by x starts to overlap the tail part of the max-PDF coincides well with the connection point e* 
for the theoretical w-PDF. 3) The center part of the min-PDFs can be adjusted, quite accurately, by the 



scaled PDFs defined by (191 with a scale factor a. It is attractive that the value of a is almost common to 
every min-PDFs with different values r/rj. 

There is no theoretical prediction yet, which is based on an ensemble theoretical aspect or on a dynamical 
aspect starting with the N-S equation, to produce the formula for the center part PDF that represents the 
contributions both of the coherent turbulent motion providing intermittency and of incoherent fluctuations 
(background flow) around the coherent motion. The discoveries given above may provide us with a correct 
pathway to formulate a dynamical theory which produces, properly, the formula for the center part of PDFs 
starting with the N-S equation. The discoveries open a new door to separate the two elements of turbulence, 
i.e., the coherent motion and the incoherent motion, and may lead us to a appropriate new method to make 
each element visible, separately, in the near future. A study to this direction is now in progress, and will be 
reported elsewhere in the near future. 

Let us close this paper by noting the following comments. It was found that the new scaling relation ^ 
is deeply related to the 5-scale Cantor sets [1] created from d°° periodic orbits [26 . In this respect, it may 
be reasonable to interpret that n introduced in ([l]) represents the number of stages in the (5-scale Cantor sets 
since n increases approximately by one for every S (see Table [I]) . On the other hand, we observe that the h 
is independent of 5, and therefore it may be appropriate to interpret that n is a good number representing 
a number of steps associate with the energy cascade model. 
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